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We present explicit expressions for the high-frequency asymptotic behavior of electron self-energy 
of general quantum impurity models, which may be useful for improving the convergence of dy- 
namical mean-field calculations and for the analytic continuation of the electron self-energy. We 
also give results, expressed in more physical terms, for the two-orbital and three-orbital rotationally 
invariant Slater-Kanamori interactions, in order to facilitate calculations of transition metal oxides. 
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Quantum Monte Carlo solutions of quantum impu- and 
rity modelsi play an important role in the dynamical 
mean-field theory (DMFT) of materialsi^ In these cal- 
culations, knowledge of the high-frequency behavior of 
the electron self-energy is useful, both to improve the 
convergence of the dynamical mean-field self-consistency 
process^"— and to provide the normalization needed for 
analytic continuation of the electron self-energyii While 
the general procedure for deriving the high frequency be- 
havior is knowUfii^rJ^ explicit expressions seem not to 
have appeared except for the single-orbital Anderson im- 
purity model^rii and its cluster generalizations J^^i^ 

In this Brief Report, we present these results for gen- 
eral quantum impurity models. We also present the 
specialization to the two and three orbital models with 
rotationally invariant (Slater-Kanamori)i''<-§ on-site in- 
teractions, relevant to late and early transition metal 
oxidesi^^— where the physics is dominated by Cg and f2g 
orbitals respectively. In supplementary material (avail- 
able online) we present a Mathematica notebook which 
generates the needed terms and presents the resulting 
output. 

A quantum impurity model involves a set of "impurity" 
degrees of freedom coupled to a non-interacting bath. Wc 
begin by defining a general impurity model. We denote 
the impurity states by an index a representing spin and 
orbital degrees of freedom and define the operator that 
creates an electron into impurity state a as Vl- We de- 
note the bath states by an energy/momentum label k and 
a spin/orbital index a, with c|.^ the operator creating an 
electron into one of these state. The impurity model is 
the sum of a non- interacting (i^nonint) and interacting 
(Hint) terms 



i/i„t= Y. /"^"^'^'^^^l^bi^b 



(3) 
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The electron self-energy E is then a matrix in space 
of impurity model states a with components E"^. An 
expression for the high-frequency behavior can be ob- 
tained from examination of the short-time behavior of 
the electron Green's function,— i^^iS which in turn can be 
expressed in terms of the commutator of the Hamiltonian 
with electron operators. Evaluating Eqs. (215) and (216) 
of Ref. Q we find: 



with 



and 



lim E'^^(z) = S;5,^ + ^, 

z— s-oo z 
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Sf= J2 ^ai'a26ifc2(^li<^biV'6.) 
aia2bifc2 

aibi A 
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where 
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In these expressions the expectation values must in gen- 
eral be evaluated from the numerical solution of the quan- 
tum impurity model. 

For the single-orbital Anderson impurity model,— the 
nonzero elements of the tensor /"i''2fcifc2 j^^^y \^q chosen 
as /tut ^ Jim ^ u/A and /titt = jim = ^u/A. One 



can readily verify that Eqs. ©-((S)) yield the well-known 
results that2r2i 
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while T,^^ can be found by flipping all spin indices, and 
all "off-diagonal" self-energies (e.g. T^^) are zero. 

While these expressions [Eqs. ©-([S])] may be straight- 
forwardly programmed and evaluated for general impu- 
rity models and interactions, it is useful to present results 
for the important special case of single-site dynamical 
mean-field theory of d electrons subject to rotationally 
invariant Slater-Kanamori interactions f^^^ because the 
results can be expressed in terms of physically relevant 
densities and the symmetries are automatically imple- 
mented, leading to fewer expressions to be evaluated. 

The on-site interacting part of the Slater-Kanamori 
Hamiltonian is traditionally expressed in terms of three 
parameters, U, U and J, as 

a a>f3,a 

+ ^ ([/' - J)n^^np, (11) 

a>fS,a 

where Uaa = i'aa'^^aa and a labels the orbitals. 

The five orbital states of the d-manifold are degener- 
ate in free space, but in cubic symmetry are split into 
a doublet, transforming as the eg representation of the 
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cubic group, and a triplet, transforming as the t2g repre- 
sentation. In many materials one may restrict attention 
to one of these two manifolds. Symmetry considerations 
then imply that the self-energy is diagonal in the orbital 
basis and, in paramagnetic states, diagonal in the spin 
states as well. 

In the case of two-orbital model (cg system), a,/3 G 
{1,2}. We can express the high-frequency behavior for 
the self-energy of a spin up electron on orbital 1 in 
terms of expectation values of orbital densities and spin- 
exchange operators as 



^'y' = U{n,0 + L/'(n2,) + {V J){n,,), 



(12) 



S^'"' = U'inuKl - inn)) + (t/')'("2i)(l - {n2i)) 
+ {U'-jf{n2t){l-{n2t)) 
+ 2UU'{{nun2i) ^ {nu){n2i)) 
+ 2U{U' - J)((ni;n2t) - {nn){n2t)) 
+ 2U'{U' - J)((n2tn2i) - {n2^){n2i)) 
+ J2 ((nu) + (n2i) - 2(ni^n2i)) 

+ 2 J( J - U) (vItV'^^V'U^2t) 

+ 2J{U + J- 2U') (^4^Vl^V'2i^2t) , (13) 

and the forms for the other orbital/spin can be straight- 
forwardly found by interchanging indices. Note that the 
off-diagonal elements are zero. 

For the three-orbital case {t2g system), a, j3 run from 1 
to 3. The high-frequency behavior of the self-energy for 
a spin up electron on orbital 1 can be similarly expressed 
as 
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: t/(nu) + [/'((n2i) + (nsi)) + iU' - J)((7i2t) + (n3t)), 

: C/2(nu)(l - (nu)) + {U'f [((n2;) + (n3^))(l - (^2;) - (ng;)) + 2(7.2^^34)] 

+ ([/' - Jf [((n2t) + (n3t))(l - {n2t) - ("3t)) + 2(n2tn3t)] + J\2{m{) + (^2^) + {n^^)) 

+ 2{UU' - J^){{nun2i) + (nunai)) - 2UU'{nu){{n2i) + (risi)) 

4- 2U{U' - J){{niin2t) + ("un-at) " ("-u)(»^2t) - (»^u)("-3t)) 

+ 2U'{U' - J) [{n2fn2l) + (n2i»^3t) + ("-2t»^3i) + ("-3t»^3i) - (("-2t) + ("-3t)) ((»T-2i) + (»^3i))] 

+ 2J(J - u) ((7At^V^^V'uV^2t) + (v4tV'li^uV'3t)) 

+ 2J{U + J- 2U') ((V'ItV'I^^2iV'2t) + (V'ItV'l4^3i^3t)) 



(14) 
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^21^2t^3iV'3t 



V'24,V'lt^2tV'3i 



^2t'^3i^2i^3t) + (^'3iV'iltV'2iV'2t 



(15) 
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while the forms for other orbitals/spins can be found by 
permuting indices. The off-diagonal elements are also 
zero. In this formulation only 27 correlators (with four 
fermionic operators) need to be measured, rather than 



the 6* = 1296 correlators arising in a straightforward 
implementation of the general equations above. 

As the number of orbitals increases, the number of 
terms proliferates, and for example in a treatment of the 



full 5-fold degenerate d manifold more than 100 terms 
occur. The expressions are too lengthy to present in the 
published version, but we include as supplementary in- 
formation a Mathematica notebook which generates the 
needed terms and presents the resulting output. In the 
same file we have also presented the calculation for the 
two and three orbital cases. This code can also easily be 
manipulated to treat Hamiltonians with arbitrary inter- 
actions. 
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